arXiv:1507.05177vl [astro-ph.EP] 18 Jul 2015 


DES 2015-0052 
Fermilab PUB-15-299-AE 


Observation of Two New L4 Neptune Trojans in the Dark Energy 

Survey Supernova Fields 


D. W. Gerdes\ R. J. Jennings^, G. M. Bernstein^, M. Sako^, F. Adams^’^, D. Goldstein^’®, 
R. Kessler^’8, T. Abbott^ F. B. Abdalla^o, S. Allamii, A. Benoit-Levy^o, E. Bertin^^-is^ 
D. Brooks^°, E. Buckley-Geer^^, D. L. Burke^^’^®, D. Gapozzi^®, A. Garnero RoselF^’^®, 
M. Garrasco Kind^®’^°, J. Garretero^^’^^, G. E. Gunha^^, G. B. D’Andrea^®, L. N. da 
Gosta^^’^®, D. L. DePoy^^, S. Desai^^’^^, J. P. Dietrich^^’^®, P. DoeP°, T. E. Eifler^’^^, 

A. Fausti Neto^^, B. Flaugher^^, J. Frieman^^’^, E. Gaztanaga^^, D. Gruen^®’^®, 

R. A. GruendP®’^®, G. Gutierrez^^, K. Honscheid^®’^®, D. J. James®, K. Kuehn^^, 

N. Kuropatkin^^, O. Lahav^®, T. S. LP^, M. A. G. Maia^^’^®, M. March^, P. MartinP®’^^, 
G. J. MilleG’\ R. Miquep2>3^ R. G. NichoF®, B. Nordl^ R. Ogando^^’is, A. A. Plazas^^ 
A. K. RomeP^, A. Roodman^^’^®, E. Sanchez®®, B. Santiago®®’^^, M. SchubnelF, 

L Sevilla-Noarbe®®’^®, R. G. Smith®, M. Soares-Santos^\ E. Sobreira^^’^^, E. Suchyta^®’®®, 
M. E. G. Swanson^®, G. Tarle^, J. Thaler®^ A. R. Walker®, W. Wester^^, Y. Zhang^ 

(The DES Gollaboration) 


ABSTRACT 

We report the discovery of the eighth and ninth known Trojans in stable orbits around Nep¬ 
tune’s leading Lagrange point, L4. The objects 2014 QO 441 and 2014 QP 441 were detected in 
data obtained during the 2013-14 and 2014-15 observing seasons by the Dark Energy Survey, 
using the Dark Energy Camera (DECam) on the 4-meter Blanco telescope at Cerro Tololo Inter- 
American Observatory. Both are in high-inclination orbits (18.8° and 19.4° respectively). With 
an eccentricity of 0.104, 2014 QO 441 has the most eccentric orbit of the eleven known stable 
Neptune Trojans. Here we describe the search procedure and investigate the objects’ long-term 


dynamical stability and physical properties. 
Subject headings: minor planets, asteroids: general 
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Trojan asteroids, together with main belt as¬ 
teroids and members of the classical Kuiper Belt, 
constitute the only dynamically stable populations 
of minor planets in the solar system. Trojans share 
the same orbital period as a major planet, leading 
or trailing the major planet by approximately 60 
degrees in an orbit centered upon the L4 or L5 
Lagrange point. Trojan asteroids are associated 


with a number of planets, including Earth (Con- 
et al.|2011| ). Jupiter has the largest and best- 


nors ( 


characterized population, with over 6000 known 
Trojans. The total number of >1 km-sized Jupiter 


Trojans is estimated to exceed 600,000 (Yoshida & 


Nakamura 2005), comparable to the similar-sized 
population of the main asteroid belt ( [Jewitt et al. 


2000). The number of large (> 65 km) Neptune 


Trojans may exceed their Jovian counterparts by 


more than an order of magnitude (Chiang & Lith- 


wick||2005 Sheppard fc Trnjillo||2006). 

Since the discovery of the first Neptune Tro¬ 


jan, 2001 QR 322 (Chiang et al. 2003), only eight 


additional Neptune co-orbitals have been discov¬ 
ered prior to this workQ Seven of the nine previ¬ 
ously known Neptune Trojans occupy the L4 re¬ 
gion. There is no a priori reason to expect the 
L4 and L5 populations to be particularly different 
( Sheppard fc Trujillo||201Q ); the dearth of known 
L5 Trojans can be ascribed to the fact that Nep¬ 
tune’s L5 region is presently on a line of sight to 
the galactic center, a crowded field of point-like 
sources against which the detection of transients 


is observationally challenging (Sheppard & Tru- 
jillo||2QlQ Parker et al.||2Q13 ). 


In addition to their intrinsic interest, the popu¬ 
lation of Neptune Trojans provides an important 
set of constraints on the dynamical evolution of 
our Solar System. Although these bodies have 
relatively little total mass, they play the role of 
a canary in the coal mine: such small bodies act 
as test particles and are easily disrupted through a 
variety of channels, especially in the early epochs 
of Solar System history. Possible mechanisms that 


^The Minor Planet Center lists a 9th object, 2004 KVis, 
as an L5 Trojan, but with an eccentricity of 0.18 and a 
libration amplitude of ~ 70°, it is known to be unstable 
on ~Myr timescales and is likely a temporarily captured 
scattered disk object ([Nesvorny fc Dones||2002 [Horner fc| 


Lykawka|2012 Guan et al.|2012|. 
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could affect the orbits of these bodies include gas 


drag in the early solar nebula (Peale 1993 Mur- 
ray||1994 ), planetary migration ( Gomes||1998 


Ko- 


rtenkamp et ah 2004), the growing mass of plan¬ 
ets during their formation (Marzari fc Scholl|1998 


Fleming fc Hamilton|2000 ), encounters with pass¬ 

ing stars in the solar birth cluster (which can easily 
perturb Neptune’s orbit, see Li & Adams 2015), 


and a wide range of resonance phenomena (Mor- 


bidelli et ah 2005). Determining the properties 


and orbital characteristics of the Trojan popula¬ 
tion also places interesting constraints on theories 
of their formation ( Chiang fc Lithwick||2005 ). For 
example, Parker (2015) uses the high mean incli¬ 


nation of Neptune Trojans to argue that the disk 
into which Neptune migrated must have been dy¬ 
namically excited prior to Neptune’s arrival. 


In this paper, we report the discovery of the 
eighth and ninth known L4 Neptune Trojans, 
2014 QO 441 and 2014 QP 441 . This paper is or¬ 
ganized as follows. Section briefly describes the 
Dark Energy Survey and camera, DECam. In Sec¬ 
tion we describe the supernova data reduction 
pipeline that is the source of transient candidates 
considered in this analysis, discuss our method 
for identifying distant solar system objects, and 
present our observations of the two Neptune Tro¬ 
jans. In Section we describe numerical investiga¬ 
tions of their long-term dynamical stability. Sec¬ 
tion describes their physical characteristics in 
the context of the broader population of Trojans. 
We conclude in Section with a summary of our 
results and a discussion of their implications. 


2. The Dark Energy Survey 


The Dark Energy Survey (DES; Elaugher|2QQ5 ) 
is a five-year optical imaging survey being car¬ 
ried out on the 4-meter Blanco telescope at Cerro 
Tololo Inter-American Observatory in Chile. Ob¬ 
servations are car ried out with the Da rk Energy 
Camera (DECam; Elaugher et al.|[2QT5 ), a 3 deg^ 
prime-focus camera whose focal plane consists of 
62 2kx4k fully-depleted, red-sensitive CCDs. DE¬ 
Cam saw first light in September 2012 , and un¬ 
derwent “Science Verification” operations between 
November 2012 and Eebruary 2013. DECam’s 
power as a discovery instrument for distant so¬ 
lar system objects was demonstrated with the dis¬ 
covery, during community time early in the com¬ 


missioning phase, of the Sedna-like dwarf planet 
2012 VP 113 ( Trujillo fc Sheppard||2014~ ), which 
has the most distant perihelion (> 80 AU) of any 
known solar system object. Eull-scale survey op¬ 
erations began in August 2013 and will continue 
through at least 2018. 

To achieve the goal of measuring the dark en¬ 
ergy equation of state to high precision, DES is 
divided into two distinct, interleaved surveys. The 
Wide Survey covers 5000 deg^ of the south galac¬ 
tic cap in the grizY bands, imaging each survey 
tile approximately ten times in each band over 
the course of the survey in order to perform high- 
statistics measurements of weak gravitational leas¬ 
ing, galaxy-galaxy correlations, and properties of 
galaxy clusters. The DES Supernova Program 
(DES-SN, Bernstein et al. ( 2012 | )) images ten dis¬ 
tinct DECam fields (8 “shallow” and 2 “deep”) in 
the griz bands at approximately weekly intervals 
throughout the DES observing season, which runs 
from mid-August through mid-Eebruary. Charac¬ 
teristics of these fields are shown in Table m Al¬ 
though primarily intended to find and characterize 
iarge numbers of Type la supernovae, the super¬ 
nova fieids are also excellent hunting grounds for 
distant minor planets, which move slowly enough 
that they can remain in the same field of view for 
weeks or months at a time, and even from one 
DES observing season to the next. In particular, 
the two Stripe-82 fields (“S-fields”) and the three 
XMM-LSS fields (“X fields”) have moderate eclip¬ 
tic latitudes ranging from —20 to —15 degrees, and 
also, serendipitously, lead the present position of 
Neptune by approximately 60 degrees of ecliptic 
longitude. We search for Neptune Trojans in these 
fields. 


3. Observations and Search Strategy 


This analysis uses imaging data processed for 
the Supernova Survey by the DES data manage¬ 
ment (DESDM) pipeline ( Mohr et al.||2012 Desa 
et al.| 2012 | ) at the National Center for Super¬ 


computing Applications (NCSA). After cleanup 
and detrending, astrometric solutions are ob¬ 
tained using SCAMP (Bertin 2006| ), with stellar 


positions in each image matched to the UCAC4 
catalog ( Zacharias et al.|| 2012 ) to an absolute as¬ 
trometric precision of roughly 0.1". Image sub¬ 
traction is performed with the Hotpants algo- 
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Table 1 

Summary of DES supernova fields. 


Field 

Depth 

Center (equatorial) 

RA (hours) DEC (deg.) 

Center (ecliptic) 

A (deg.) /3 (deg.) 

Cl 

shallow 

03:37:05.83 

-27:06:41.8 

42:49:30.2 

-44:52:29.6 

C2 

shallow 

03:37:05.83 

-29:05:18.2 

41:52:47.4 

-46:44:19.2 

C3 

deep 

03:30:35.62 

-28:06:00.0 

40:25:56.8 

-45:19:21.4 

XI 

shallow 

02:17:54.17 

-04:55:46.2 

30:28:17.9 

-17:38:53.3 

X2 

shallow 

02:22:39.48 

-06:24:43.6 

31:06:23.4 

-19:26:48.6 

X3 

deep 

02:25:48.00 

-04:36:00.0 

32:31:58.0 

-18:00:28.3 

SI 

shallow 

02:51:16.80 

00:00:00.0 

40:22:15.4 

-15:41:10.2 

S2 

shallow 

02:44:46.66 

-00:59:18.2 

38:26:48.0 

-16:07:37.7 

El 

shallow 

00:31:29.86 

-43:00:34.6 

346:05:17.0 

-41:44:05.3 

E2 

shallow 

00:38:00.00 

-43:59:52.8 

346:43:33.9 

-43:11:58.4 


rithm (Becker 2015) using PSF-matched deep 
templates obtained during Science Verification. 
Source detection is performed on the subtracted 


images using Sextractor (Bertin & Arnouts 


1996). A machine-learning algorithm (Goldstein 
et al.||2Q15 ) reduces artifacts from instrumental 


or 


image-subtraction effects by a factor of nine. All 
data have been reprocessed using the difference¬ 


imaging pipeline described in Kessler et al. (forth¬ 


coming). The efficiency of the pipeline for recon¬ 


structing and accepting true point-like sources is 
measured by inserting fake supernovae into the 
images using the SNANA code ( [Kessler et al.|2QQ9 ) 
prior to image-subtraction. The efficiency reaches 
50% for signal-to-noise ratio (SNR) of about 5 (see 
Fig 7 in Kessler et al. ( forthcomingj )), which corre¬ 
sponds to an average magnitude depth of 23.5 in 
the shallow fields and 24.5 in the deep fields. We 
compute PSF magnitudes in the AB system. Im¬ 
ages are calibrated nightly to tertiary standards 
to a photometric precision of 2 % on average. 

The r- and i—bands have the best S/N, with 5 cr 
limiting magnitudes of 23.8 in the shallow fields 
and 24.5 in the deep fields. We therefore carry 
out the initial candidate search in these bands, 
and add observations from g and z-band once can¬ 
didates have been identified. Our target fields 
contain approximately 185,000 single-epoch tran¬ 
sients in the ri bands from the first (Aug. 2013 to 
Feb. 2014) and second (Aug. 2014 to Feb. 2015) 
DES observing seasons. 

An object at 30 AU undergoing retrograde mo¬ 
tion moves at an apparent rate of up to 4" per 
hour. In the shallow fields, exposure times range 
from 150 s in r-band to 400 s in z-band. These 
exposures are sufficiently short that such an ob¬ 


ject appears stellar and is accepted with high ef¬ 
ficiency by the machine-learning cut. In the deep 
fields, however, image subtraction and source de¬ 
tection are ordinarily performed on co-added se¬ 
quences that range from 3 x 200 s in ^-band to 
11 X 360 s in 2 ;-band. Because such a procedure 
would frequently result in the rejection of slow- 
moving transients, for the deep fields we have run a 
modified version of the transient detection pipeline 
that uses the single-epoch exposures. This pro¬ 
cedure has the added advantage of producing a 
fine-grained time series well-suited for light-curve 
measurements as discussed in Section [H 


The apparent motion of a distant solar system 
object in the roughly one week between DES visits 
to each field is due largely to Earth’s reflex motion. 
Beginning with the list of transients identified in a 
given visit by the difference-imaging pipeline, we 
search for counterparts at two subsequent visits in 
a window consistent with seasonally-appropriate 
reflex motion at a rate of less than 150"/day. 
Triplets selected via this geometrical technique are 
tested for goodness of fit to an orbit using code 
derived from the work of IBernstein fc Khushalaml 
( 2000 ). Observations from other bands and visits 
are added iteratively. 

2014 QO 441 and 2014 QP 441 were initially de¬ 
tected at apparent r-band magnitudes of 23.3 and 
23.9 in the deep XMM-LSS field on Aug. 21 , 2014, 
at ecliptic latitudes of —18.5° and —18.0° respec¬ 
tively. 2014 QO 441 was observed 69 times on 
12 different nights through Jan. 29, 2015, while 
2014 QP 44 iwas observed 62 times on 13 nights 
through Jan. 22 , 2015. In addition, 2014 QO 441 
was recovered in five 90-second wide-survey expo¬ 
sures from Oct. 2 and 13, 2013. The trajectories of 
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these two objects are shown in Figure Orbital 
elements and other properties of these two objects 
are summarized in Table |2j The 2013 observations 
of 2014 QO 441 increase its arc length by nearly a 
year and substantially reduce the uncertainty in 
the orbit. A summary of these objects’ proper¬ 
ties in comparison with the other known Neptune 
Trojans is shown in Table 


To compute the libration amplitudes and peri¬ 


ods we follow the procedure of Parker et ah (2013) 


and Parker (2015). Using the same 1000 clones of 


each object, we define the half-peak rms libration 
amplitude for each clone to be 


N 


1/2 


^11= I 


( 1 ) 


i=l 



Fig. 1 .— Observed trajectories of 2014 QO 441 
(lower curve) and 2014 QP 441 (upper curve) rel¬ 
ative to the three DES XMM-LSS supernova 
fields over the first two DES observing seasons. 
The cyan portions of the trajectories encompass 
the 2014-15 DES observing season. Large dots 
indicate nights on which the objects were de¬ 
tected. Each rectangle represents one DECam 
CCD. 2014 QO 441 was also detected in wide- 
survey exposures outside of the supernova fields 
on two nights in October, 2013. 


4. Dynamical Stability 


To investigate the long-term behavior of the 
orbits, we have simulated swarms of 1000 clones 
each of 2014 QO 441 and 2014 QP 441 . The initial 
orbital elements of each clone were drawn from a 
multivariate normal distribution using the full six¬ 
dimensional covariance matrix obtained from the 
fit to the observations. We followed each clone for 
10 Myr in the presence of the four giant planets 
using the hybrid symplectic algorithm in the Mer- 
CURy6 Wbody integrator ( Chambers|1999 ). Each 
of the clones remained in a bound orbit about L4, 
evolving stably in aei-space. An example is shown 
in Eigurej^ 


where i runs over the N samples of each clone 
during the 10 Myr integration. The resonant ar¬ 
gument 011 = Atv — At, where A at and At are the 
mean longitude of Neptune and the Trojan respec¬ 
tively (Chiang et al. 2003). The mean longitude 
computed at each epoch of the integration is 

A = M + U + u;, ( 2 ) 


where M is the mean anomaly, Q is the longitude 
of the ascending node, and uj is the argument of 
perihelion. Libration periods are obtained from a 
spectral decomposition of 0 (t), and errors are ob¬ 
tained from the rms scatter among the clones. We 
find that 2014 QO 441 librates with an rms ampli¬ 
tude of 10.9 ± 0.1° and a period of 9074 ± 3 years, 
while for 2014 QP 441 we find a libration amplitude 
of 7.8 ± 1.3° and a period of 9114 ± 5 years. The 
orbital motion of both Trojans in the co-rotating 
Neptunian frame over the course of one full libra- 
tion period is shown in Eigurej^ 

In addition to the ^9100-year libration peri¬ 
ods, we observe oscillations in eccentricity and in¬ 
clination on 1-2 Myr timescales, which are a by¬ 
product of the long-term exchange of angular mo¬ 
mentum between Jupiter and Saturn. To inves¬ 
tigate these further we have continued 32 of the 
integrations for each object for 1 Gyr. We observe 
that these oscillations continue, but do not grow 
in amplitude, over the full timespan of the inte¬ 
grations for 2014 QP 441 . Eive of the 32 clones of 
2014 QO 441 did not survive the full 1 Gyr inte¬ 
gration, with the earliest instability occurring be¬ 
tween 500-600 Myr. Thus 2014 QP 441 appears to 
be completely stable in its present configuration, 
while 2014 QO 441 possesses a dynamical half-life 
on the order of 4 Gyr. Both objects could there¬ 
fore be primordial members of the disk into which 
Neptune migrated. 
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Table 2 


Orbital elements and other properties of the two new Trojans. 


Parameter 

2014 QO 441 

2014 QP 441 

a (AU) 

30.102 ± 0.001 

30.108 ± 0.011 

e 

0.1046 ± 0.0004 

0.067 ± 0.002 

i (deg) 

18.8316 ± 0.0003 

19.405 ± 0.001 

u (deg) 

113.1 ± 1.0 

1.74 ± 0.67 

Q (deg) 

106.998 ±0.004 

96.538 ± 0.007 

Perihelion date 

1937/11/28 ± 176d 

2043/4/16 ± 79d 

Epoch JD 

2457000.5 

2457000.5 

Libration 
period (yr) 
Libration 

amplitude (deg.) 

9074 ± 3 

9114 ± 5 

10.9 ± 0.1 

7.8 ± 1.3 

Arc length (days) 

484 

154 

Apparent mag (r) 

23.3 

23.9 

Hv 

8.2 

9.1 

r — i (mag.) 

0.21 ± 0.10 

0.20 ± 0.06 

i - z (mag.) 

0.01 ± 0.11 

0.08 ± 0.08 

Diameter (km)®^ 

120 

70 


Assuming 5% albedo. 


Table 3 

Properties of Neptune Trojans 


Name 

a (AU)^ 

e 

i (deg.) 

Hv 

L 

diam. (km)’^ 

Libration ampl. (deg.)*^ 

Libration per. (yr.) 

2001 QR 322 

30.207 

0.027 

1.3 

8.2 

L4 

140 

25.5tli 

9200 ± 3 

2004 UPio 

30.130 

0.024 

1.4 

8.8 

L4 

100 


8874 ± 8 

2005 TN 53 

30.127 

0.067 

25.0 

9.0 

L4 

80 


9428 ± 6 

2005 TO 74 

30.120 

0.055 

5.3 

8.5 

L4 

100 

9 2”^°’^ 

^•^- 0.5 

8822 ± 10 

2006 RJios 

30.045 

0.032 

8.2 

7.5 

L4 

180 

6 3+°-^ 

O-'J_0.3 

8858 ± 10 

2007 VL 305 

30.079 

0.065 

28.1 

8.0 

L4 

160 

-‘-^•^- 0.10 

9619 ± 8 

2008 LCis 

29.922 

0.085 

27.6 

8.4 

L5 

100 


9516 ± 7 

2011 HM 102 

30.059 

0.080 

29.4 

8.1 

L5 

140 

P 

bo 

1 + 
0 0 

9545 ± 4^ 

2012 UV 177 

30.032 

0.074 

20.8 

9.2 

L4 

80 

13® 

9334 ± 10 

2014 QO 441 

30.089 

0.105 

18.8 

8.2 

L4 

120 

d d 
+ 1 
P 
d 

9074 ± 3 

2014 QP 441 

30.108 

0.067 

19.4 

9.1 

L4 

80 

7 9“^^'^ 

‘ — 1.3 

9114 ± 5 


®^Barycentric osculating elements at epoch 2014/12/09. 
'^Assuming 5% albedo. 

'^Values for first 8 objects from |Parker| |20T^ . 

" ^Parker et al.| ( |2013| ) 

^Alexandersen et al.| ( |2014| ) 
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Fig. 2.— Evolution of the resonant angle and 
the orbital parameters e and i over a lOMyr inte¬ 
gration for the best-fit orbits of 2014 QO 441 (top) 
and 2014 QP 441 (bottom). This pattern continues 
stably for each of the 1000 clones analyzed. 



2014 Q0441 
2014 QP441 


10 15 

X (AU) 


Fig. 3.— Orbital motion of 2014 QO 441 (blue) 
and 2014 QP 441 (red) in the ecliptic plane over 
the course of one full ^ 9100-year libration period 
as seen from a frame rotating counterclockwise at 
the speed of Neptune. In this coordinate system 
the Sun is at (0,0) and Neptune lies along the x- 
axis between x = 28.9107 and x = 30.3271. Each 
individual sub-cycle of the full pattern corresponds 
to one 165-year orbit. 


5. Physical Characterization 


The measured r — i and i — z colors shown in 
Table [^indicate that both objects are slightly red. 
The derived V — R magnitudes of ^ 0.44±0.11 are 
consistent with previously reported colors of the 
L4 and L5 Neptune Trojans ( Parker et al.|2013 ) as 
well as Jupiter Trojans and neutral Centaurs, and 
are slightly bluer than classical KBOs. This sim¬ 
ilarity may indicate a common formation mech¬ 
anism and history, although future near infrared 
photometry may reveal otherwise unapparent pop¬ 
ulation differences as is known to be the case with 
Jupiter Trojans ( Wong et al.||2014 ). 

We examined the light curves of each object 
to search for periodic behavior indicative of ro¬ 
tation or a possible binary system. We carried 
out this analysis using both the Phase Dispersion 


Minimization (Stellingwerf 1978) and generalized 


Lomb-Scargle periodogram (Scargle 1982 Zech- 


meister fc Kurster|2QQ9 ) techniques. We found 


no 


evidence for modulations on timescales between 
0.5-48 hours, leading us to conclude that within 
our photometric sensitivity both 2014 QO 441 and 
2014 QP 441 are relatively round and featureless. 
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6. Summary 


This paper reports the discovery of 2014 QO 441 
and 2014 QP 441 which represent the eighth and 
ninth members of the L4 family of Neptune Tro¬ 
jans. They appear to be dynamically stable on 
> Gyr timescales and have colors similar to other 
Trojans. Their high inclinations lend further sup¬ 
port to the idea that Neptune Trojans are a dy¬ 
namically “hot” population that may have been 


excited prior to capture by Neptune (Parker 2015) 


The Dark Energy Survey will have further oppor¬ 
tunities to observe 2014 QO 441 during the 2015- 
16 and 2016-17 seasons, when it will again enter 
the deep XMM-LSS supernova field. 2014 QP 441 
will also be visible in this field during the 2015-16 
season, then will briefly pass through one of the 
Stripe-82 fields during the 2016-17 and 2017-18 
seasons. Further measurements of this object will 
improve the estimate of its dynamical lifetime. 

The existence of the Neptune Trojans, along 
with other small bodies in the Solar System, poses 
a set of interesting dynamical problems. These 
bodies could have been placed in their current or¬ 
bits at the epoch of solar system formation, or per¬ 
haps more recently. In the former scenario, the 
orbits must be stable over the ^ 4.6 Gyr age of 
the Sun. Sufficiently large perturbations to the 
orbit of Neptune, or the other giant planets, could 
lead to the removal of Trojans from their libra- 
tions about the Lagrange points (e.g.. Figure 
shows that small perturbations of Neptune due to 
the other planets lead to changes in the resonant 
angle). As a result, the long-term stability of these 
orbits implies upper limits on the degree of disrup¬ 
tion suffered by the solar system. On the other 
hand, if the placement of the Trojans into their 
currently observed orbits is more recent, then the 
mechanism that populates such orbits must be un¬ 
derstood. These issues thus present a rich set of 
open dynamical questions. 

Although this work has emphasized the detailed 
study of these two new Trojans, an additional 20 
Kuiper Belt and other trans-Neptunian objects 
discovered in our search of the DES supernova 
fields have been given provisional designations by 
the Minor Planet Genter. The most notable of 
these is 2013 RFgg, which has a semi-major axis 
of 325 AU and a perihelion of 36 AU, making it 
one of the longest-period trans-Neptunians known. 


These objects will be described in forthcoming 
work. 

With three years remaining in the Dark En¬ 
ergy Survey, we can confidently anticipate that 
several dozen new trans-Neptunian objects, and 
perhaps including more Neptune Trojans, will en¬ 
ter the supernova fields and be discovered by the 
techniques described in this paper. But the su¬ 
pernova fields represent less than 1 % of the full 
5000 sq. deg. wide survey area. Although the 
sparse cadence with which a given wide survey 
region is observed—typically 2-4 times per sea¬ 
son in each of five filters—makes the identifica¬ 
tion of slow-moving transients more challenging, 
a search in the wide survey fields is likely to re¬ 
sult in hundreds of TNOs to a limiting magni¬ 
tude of r ^ 23.8. Moreover, the wide survey foot¬ 
print includes a small region near the ecliptic, but 
is distributed primarily over higher ecliptic lati¬ 
tudes, an area not well covered by previous TNG 
searches. This makes DES particularly well suited 
for the detection of the higher-inclination, “hot” 
TNG populations. DBS’s combination of area and 
depth suggests that DES can discover 6-10 times 
more such objects than any previous survey. Be¬ 
cause the detached / inner Gort cloud population 
is even less likely to be confined to the ecliptic re¬ 
gion, we expect DES to have at least an order of 
magnitude more sensitivity to this very interest¬ 
ing and rare group of objects compared to prior 
searches, shedding further light on the processes 
that shaped the solar system. 
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